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Abstract 



The Born-Infeld form of the hydrogen atom has a spectrum that can be used to determine the 
physical viability of the theory, and place an experimentally relevant bound on the single parameter 
found in it. We compute this spectrum using the relativistic Dirac equation, and a form of the Born- 
Infeld potential that approximates the self-field corrections of the electron. Using these together, we can 
establish that if the Born-Infeld nonlinear electrodynamics is to be physically relevant, it must contain 
a fundamental constant that is well-below the original value proposed by Born. This work extends 
the original Schrodinger spectrum from [T] for the self-field correction, and shows that using the Dirac 
equation introduces minor corrections - but also gives access to a range for the fundamental constant 
that is below that attainable from non-relativistic considerations. 



1 Introduction 

The motivations for studying Born-Infeld electrodynamics [3] are varied - while the original motivation 
(finite point-source self-energy) may or may not be compelling in modern times, the utility of the theory as, 
at the very least, a proving ground for potential nonlinear modification is not to be undervalued Of 
particular interest here is the role of particles, which both source and respond to nonlinear fields [8]. In 
addition to providing a vehicle for thinking carefully about nonlinearity (in preparation for other nonlinear 
studies, like GR, say [5]), the particulars of Born-Infeld have current proponents from string theory (by 
now, this is established enough to warrant a chapter in |10p. general relativity, atomic physics [111 I12| and 
spectroscopy [13] . Our work serves both as an attempt to collect and compare known spectra, as well as 
present the results of the Dirac spectrum, all in a unified language. The ultimate goal of any such spectral 
study is to place a bound on the parameter governing the strength of the Born-Infeld potential. This single 
number, the only input in the theory, can be used to evaluate the physical relevance of the Born-Infeld 
approach, in addition to setting the range of applicability of the theory. 

The Born-Infeld field equation for the electrostatic potential is: 



where a is the Bohr radius of hydrogen, a is a parameter that sets the value of the electrostatic energy 
density at r = and e is the charge of the electron. Born proposed a value of clb = g 5(1/4, 1/4) a 2 (where 
a = 7.2973525376 x 10~ 3 , the fine structure constant, and B is Euler's /3 function) that comes from equating 
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the (now finite) field energy with the rest energy of the electron (the actual number is usually given with 
dimensions of length, as — «b a). If we set the density on the right to p — e <5 3 (r), corresponding to a point 
particle, we recover the usual expression for the potential energy of a charge — e moving in the field of a 
central charge e: 



Previous work on spectra has focused on this "test-particle" approximation (as in [TH [T3]), where 
only the single-particle Born-Infeld point potential is used. This approach is justified in the limit that the 
proton charge is much larger than the electron's, so that the electron's contribution to the total Born-Infeld 
field is negligible. Of course, for hydrogen, this approximation is invalid, but the test particle potential 
is still useful. The Born-Infeld field equations are nonlinear, and so the dipole field (of a proton and an 
electron) is not the sum of the individual point fields. It is this dipole field, associated with the density 
p = e S 3 (r — r + ) — e 8 3 (r — r_), that one would ultimately like to use as a potential. In three dimensions, the 
dipole field has no closed form, and so the test particle potential is relevant to spectral studies as the only 
known exact particle solution. 

An approximation to the dipole potential was proposed [T5J that includes some of the self-field effects 
of the electron. In [T], the non-relativistic Schrodinger spectrum was calculated with this potential, and 
it was shown that in this approximation, the parameter a in the Born-Infeld equation must be less than 
Born's original value. In this paper, we extend the spectrum to include relativistic effects by computing 
the spectrum of hydrogen, with this modified potential, using the Dirac equation. By including relativistic 
effects, we can probe values of the parameter a below the original as given by Born - when we compare the 
spectra obtained by the Schrodinger and Dirac equations, the differences between the test-particle potential 
and the self-field potential appear at roughly an order of magnitude larger than the relativistic corrections 
for a = as- Using the Dirac equation allows us to distinguish between spectral differences between the two 
potentials, and those arising from relativistic effects for values of a well below as- In addition to corrections 
to the energies themselves, the Dirac spectrum will provide angular degeneracy information that can be used 
to further constrain the allowed values of a. 

We will show that as far as the test-particle form is concerned (relevant if we had a system in which 
the "electron" had charge much less than the "proton"), deviation of spectra (in all cases) from accepted 
values is minimal for "reasonable" values of the parameter found in the Born-Infeld potential (i.e. near 
the one proposed by Born). If one takes into account the electron self-field, then parametric values are 
constrained by angular degeneracy (this was shown for the non-relativistic case in [T]). After setting up the 
physical framework (the potentials), we discuss the numerical method used to find all spectra in this paper, 
the simplest method possible. Then we set the two parameters of the method, and indicate its limitations, 
by matching the Dirac hydrogen spectrum (using the Coulomb potential). From there, we are ready to 
calculate all flavors of spectra, and finally, show that in the self-field case, the angular portion of the Dirac 
spectra would be measurably different if BI holds, even for values of the BI parameter down to one-tenth of 
the theoretical limit. It is important to note that this self-field potential is still approximate - in order to 
treat the full hydrogen problem, one would need a three-dimensional dipole solution to the Born-Infeld field 
equations, absent a fully quantum form for the field. 



Here, we briefly review the self-field potential found in |15j . and set the form and dimensionless variables 
appropriate for both Schrodinger and Dirac investigations. 

The original Born-Infeld potential energy for a particle of charge e interacting with a test particle of 
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charge — e was given in Kiessling has argued [5] that when considering a two-body system like hydrogen, 
we must include the self-field of the electron, in addition to the central charge, and gives an augmented 
potential, approximating this correction, of the form: 
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We have used the identifiers 1 and 2 in ^ and Q to distinguish between these two potentials, and remind 
us of the number of charges being approximated. 

For Schrodinger's equation, we start from the separated radial ODE, with u(r) — r R(r) a function of 
r only. If we set r — pa for Bohr radius a and dimensionless p, then the radial equation can be written: 
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with u" = ^-r. The energy on the right, for the Coulomb potential, would be 
variables, the (test-particle) Born-Infeld potential for use in Schrodinger's equation is: 

dx 
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/ith similar modification for the potential ([3]). 

Using the same dimensionless p, the Dirac equation can be brought to the form: 
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where u and v make up the radial portion of the full spinor. We define k = j + 4 for j the total angular 
momentum. 



3 Method 



We use the same numerical method to find the spectra in all cases - finite differences on a regular grid 
with pj — j Ap (for some choice of step Ap). A centered difference replaces the second derivative u" in the 
Schrodinger equation, and a similar centered difference is applied to the first derivative(s) appearing in the 
Dirac equation. 

In particular, the second derivative u" appearing in Schrodinger's equation can be written approxi- 
mately as 

u"( Pj ) « + + 0(Ap 2 ) (7) 

so that the second derivative of u at the grid point pj consists of a combination of the values of the function 
u at Pj+i, Pj and P j-\. While it is clear how to embed this information in a matrix for most grid-points, the 
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"boundary points" (like pi) are more subtle. Our grid starts at pi = Ap, suppose it ends at p^ = N Ap for 
some integer N. We will need to refer to the values of the function u at zero and (N + 1) Ap to construct 
the approximation to the second derivative of u at p\ and pjv- At the left-hand boundary, the situation 
is straightforward - the u in Q is r R(r) for R(r) the radial portion of the wavefunction. For hydrogenic 
R n i(r), we know that R n e(0) is zero or finite, so w(0) = automatically. That means we can ignore the 
reference to uq appearing in the finite difference, since it is zero. 

The other end is more difficult - we expect all u to go to zero as p gets "large" (compared to unity) , so 
we choose an artificial infinity, called p^ , and demand that the radial portion of the wavefunction vanish at 
this necessarily finite value. This approximation is not so bad for low-lying energy states, where we expect 
decaying exponentials in p to dominate the behavior of u for large values of p. Still, we must be aware that 
higher energy states will have greater error associated with our choice. For now, we include our choice of 
Poo and N, the number of grid points, as inputs in our method, and then define Ap = ^ff- as the step size. 

These approximations render the ODEs as algebraic eigenvalue problems of the form: 

Bu = Eu, (8) 

for the vector u with entries Uj « u(Pj)- Our goal, then, is to solve for the energies E that are just the 
eigenvalues of the discretized differential operator D (a sparse matrix). Since we are interested in the low-end 
of the spectrum, we only want a few of the smallest eigenvalues, and we use the Lanczos iterative method 
to find just the first few eigenvalues (quickly) [j] Note that an identical procedure can be applied to the first 
order Dirac pair, although we have to be careful to include both the u and v contributions correctly - this 
means that u appearing in Q is now a vector that includes all discrete values of u and v (a vector, then, of 
length 2N), and D G R 2Nx ™. 



3.1 Test and Parameters 



As a check of the method, and a validation of the parameters poo and N used in the remainder of this paper, 
we reproduce the Dirac hydrogen spectrum, probing (successfully) both its principle quantum number and 
angular dependence by comparing with the known spectrum (using our dimcnsionlcss variables): 
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where, again, n = j + |, and n is the principle quantum number (for the Schrodinger spectrum). 



Using poo = 100, and N = 20000, we obtain a ground state energy for the Dirac equation (with the 
Coulomb potential: V(p) = — 1/p) of: E s» —1.000013313195 - this compares well with the exact value 
from ([9| - the energy values for the first three principle quantum numbers, and first two orbital states are 
shown in Table 1. Notice, there, that the lowest energy levels (n = 1, K = 1 and n = 2, k — 2) in each 
series are the most accurate - our direct method loses accuracy as energy increases, which is not surprising. 
For the usual radial wavef unctions, higher energy implies longer spatial extent, and our designation of poo 
becomes relevant in the error. In addition, the numerical spectrum is necessarily finite, so that we expect 
errors to accumulate near the "top" of the spectrum, where we ultimately approximate the infinite-energy 
scattering states with a finite value. 



practice, we could use almost any language or numerical package to implement this method. For simplicity, and 
accessibility, we used the built-in Lanczos routines found in Mathematica. 
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Table 1: Dirac spectrum calculated numerically and from ([9]). We show the numerical result up to disagree- 
ment with the theoretical values. 
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Table 2: Comparison of the spectrum computed with the test-particle potential ([2j) and the potential with 
electron self- field included ^ for Schrodinger's equation using a = as- 



4 Schrodinger Spectra 

We computed the energy spectrum as a function of a — Q clb for Q = 1 — > 30000 in steps of 1000 for 
the ground state, and first two excited states with I = and I — 1 for both the test-particle form of the 
potential ([2]), and the approximation ([3]). The results are shown on the top in Figure [IJ Note that we 
recover the two values for a that give good agreement for the hydrogen ground state energy (a hallmark of 
the potential (|3|), but it is also clear that every energy level will have two such values: At a = 0, we know 
we recover Coulomb, and we expect that for large a, the test particle potential and ^ coincide - but all the 
test-particle deviations lie above the Coulomb value, while all deviations coming from the modified potential 
start off below the Coulomb value. In order for the two spectra to match for large a, the self-field spectrum 
must come back up to the Coulomb value and cross it. 

In all cases, the potential with electron field included diverges from the Coulomb spectrum more than 
the test particle potential. In addition, the splitting of the angular degeneracy is significant when using ^ 
- it is this splitting that we will probe for values of a below Born's. In Table 2, we see the numerical values 
associated with clb for each energy level and potential. 



5 Dirac Spectra 

The general pattern for both potentials when we move to the relativistic Dirac equation is the same - very 
little changes when compared to the Schrodinger case. On the bottom in Figure [T] we see the spectra for 
K = 1, 2 — in either Born-Infeld case, the degeneracy of I = 0, 1 is split for some value of a, and the splitting 
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Table 3: Comparison of the spectrum computed with the test-particle potential ([2]) and the potential with 
electron self- field included pi using the Dirac equation (at a — as). 



is more dramatic when using ([3]). In that case, we can see an interesting crossing in the k = 2 case, where 
degeneracy is restored at some (large) value of a (near awl). 

A table of numerical values at Born's as is shown in Table 3. The data shown there is carried to the 
first digit that disagrees with the Schrodinger results in Table 2 for the potential ([3]), and enough digits to 
distinguish between the I = and 1=1 cases. 



6 Bounds on Born's Parameter 



Using the self-field potential, we see that the spectrum of hydrogen, even at Born's modest value for a, is off, 
especially when we think of the angular degeneracy associated with £ = and 1. We probe below the Born 
value, to see what type of behavior ever-smaller parametric values induce. Note that in the test particle 
case, all is well, so we focus on the potential pi). The result, for the Dirac equation, is shown in Figure [2] 



7 Conclusion 



We have demonstrated that the use of either the test-particle or modified self-field BI potential in the 
relativistic Dirac setting does not significantly change the Schrodinger results; we pick up corrections of the 
usual fine structure size in each case. Using the modified potential we see a similar divergence in angular 
dependence for the hydrogen spectra for both the Schrodinger and Dirac equations. But the difference 
between the test-particle form and the self-field form at og is approximately an order of magnitude larger 
than the relativistic corrections, hence to probe values of the parameter a below clb, and distinguish between 
the two potentials, we must use a relativistic approach. One could interpolate/extrapolate from our energies, 
as functions of a, to put concrete numerical bounds on this parameter given some experimental tolerance, 
or use our results to verify or reject an externally proposed value of a. Our relativistic calculation competes 
with the effects of spin-orbit coupling in magnitude, and it would be interesting to include some relevant 
Born-Infeld approximation to spin-orbit coupling in the current framework. 

Our spectral results exhibit I dependence that is numerically different than previous workj^J and imply 

2 In QQ, the energy for the n = 2, i = and 1=1 states are given as (in our current dimensionless scheme): E20 = .25004 
and E21 = .38202, leading to dramatic degeneracy splitting. We have checked our results for a variety of parameteric settings 
of poo and N, in addition to modifying the (necessarily numerical) integration of J3j), and in all cases (including a completely 
separate determination of the spectrum using shooting), we find agreement with the stated values in Table 2, values which are 
correctly echoed in the Dirac spectrum at a = dg as shown in Table 3. 



6 



a much milder set of spectral shifts when we use Vgj - the self-field form of the BI potential still exhibits 
greater splitting of angular degeneracy in both the Schrodinger (not shown) and Dirac settings, and this can 
be used to put a bound on the value of as as needed (depending on the application, in other words). The 
method we use is straightforward and generalizable to other nonlinear modifications, or physically inspired 
perturbations (such as spin-orbit coupling). Future work can build off of this foundation, modifying either 
Poo or N to probe different portions of the spectrum with the same accuracy. 

The authors thank D. Griffiths and M. Kicssling for valuable commentary on an early draft. 
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Figure 1: Top: the energy spectra for i — and i = 1 as computed using the test particle Born-Infeld 
potential ([2]), and the self-field potential Solid lines indicate the Coulomb-based hydrogen spectrum. 
On the bottom is the Dirac version of the same spectra, for k = 1 and k = 2, again computed using the 
test particle Born-Infeld potential, and self-field potential. The splitting of the £ = and I = 1 degeneracy 
is evident here (for example, on the bottom right, k — 2 case, the n = 3 energies begin at a — taking 
on their Coulomb value, but as a is increased, the pair becomes split in both the test-particle and self-field 
potentials). 



8 



o.nimoi (i.nrifKi: 0.00003 0.00004 o.ihhhh (i.onoof, 



-0.25nnn - 
-0.25001 '- 
-0.25002 - 
-0.25003 - 
-0.25004 - 
-0.25005 - 
-0.25006 '- 



Figure 2: The E 2 energy calculated with £ — and £ — 1 for values of a starting at a^/SO. Here, we use 
the self-field potential in the Dirac equation - the divergence of the degeneracy, for small values of param- 
eter, indicates that even for values well below Born's, the self-field Born-Infeld potential is experimentally 
unacceptable. 
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